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Abstract 

A system of periodically coupled nonlinear phase oscillators — submitted to 
both additive and multiplicative white noises — has been recently shown to 
exhibit ratchetlike transport, negative zero-bias conductance, and anomalous 
hysteresis. These features stem from the asymmetry of the stationary prob- 
ability distribution function, arising through a noise-induced nonequilibrium 
phase transition which is reentrant as a function of the multiplicative noise 
intensity. Using an explicit mean-field approximation we analyze the effect of 
the multiplicative noises being coloured, finding a contraction of the ordered 
phase (and a reentrance as a function of the coupling) on one hand, and a 
shift of the transition from anomalous to normal hysteresis inside this phase 
on the other. 
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I. INTRODUCTION 



One of the distinguishing features of the last decade has been a spectacular advance- 
ment of knowledge in the field of nanotechnology . Within a subfield that we might call 
"nanomechanics" , the topic of noise-induced transport [[I J — which concerns the mechanisms 
whereby one can extract useful work out of (nonequilibrium) fluctuations — has increasingly 
captivated researchers. 

In the early works it seemed to be a requisite for the operation of microscopic recti- 
fying devices (usually known as "molecular motors", "Brownian motors", or "ratchets") 
that — besides a built-in ratchetlike bias — the fluctuations be correlated |2|]. That require- 
ment was relaxed when "pulsating" ratchets, in which the rectifying effect arises from the 
random switching between uncorrelated noise sources, were discovered [|3|. A recent twist 
has been to relax also the requirement of a built-in bias: a system of periodically coupled 
nonlinear phase oscillators in a symmetric "pulsating" environment has been shown to un- 
dergo a noise-induced nonequilibrium phase transition, wherein the spontaneous symmetry 
breakdown of the stationary probability distribution function (SPDF) gives rise to an effec- 
tive ratchetlike potential. Some of the striking consequences of this fact are the appearance 
of negative (absolute) zero-bias conductance in the disordered phase, but near the phase- 
transition line (for small values of the bias force F, the particle current (X) opposes F), 
and anomalous hysteresis in the strong-coupling region of the ordered phase (the (X) vs F 
cycle runs clockwise, as opposed for instance to the B vs H cycle of a ferromagnet) [||]. 

Exploiting our previous experience U , in Ref . || we addressed the model using an explicit 
mean-field approach (see e.g. Ref. 0). We undertook a thorough exploration of the ordered 
phase — including the characterization of its subregions and the transition from anomalous to 
normal hysteresis in the behaviour of (X) as a function of F — and found a close relationship 
between the shape of the SPDF and the number of "homogeneous" mean-field solutions on 
one hand, and the character of the hysteresis loop on the other. 

As discussed in Ref. 0, the multiplicative noises are expected to exhibit some degree of 
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time-correlation or "colour" . Hence in this work, and as a natural continuation to Ref. || , 
we study (in mean-field approximation) the consequences of a finite correlation time r of the 
multiplicative noises in the model of Refs. . In the following sections we briefly describe 
the model and its mean- field treatment, and discuss our numerical results regarding the 
influence of colour on the phase diagram and on the transition from anomalous to normal 
hysteresis inside the ordered phase. 

II. THE MODEL AND ITS MEAN-FIELD ANALYSIS 

The model: Let us consider the following set of N globally coupled stochastic equations 
of motion, in the overdamped regime: 

dU- 1 N 

i j=i 

The stochastic variables Xi(t) are phaselike [-L/2 < Xi(t) < L/2, where L is the period 
of the Ui] and the equation is meant to be interpreted in the sense of Stratonovich. The 
second term in Eq. (|I|) models, as usual, the effect of thermal fluctuations [the are 
additive Gaussian white noises with zero mean and variance one: = , (£i(t)£j(~f)) = 

5ij5(t — £'), and T represents the temperature of the environment]. 

The "pulsating" potentials Ui(x,t) are one of the key ingredients in the model ||. In- 
cluding a "load force" F as a tool for the analysis of the noise-induced ratchet effect, their 
form is 

Ui(x,t) = V(x)+W{x)y/2Qrk{t)-Fx, (2) 

namely they consist of a static part V(x) and a fluctuating one: Gaussian noises rji(t) with 
zero mean couple multiplicatively (with intensity Q) through a function W(Xi). Whereas in 
Refs. the r)i(t) were taken as white with variance one, we now regard them instead as 
Ornstein-Uhlenbeck, with self-correlation time r: (Vi{t)Vj{t')) = ^i( <j2 /2 , r) exp(— \t — t'\/r). 

Besides being periodic, V(x) and W(x) are assumed to be symmetric V(—x) = V(x) 
and W(— x) = W(x) (there is no built-in ratchet effect). In Refs. (H| it is V(x) = W(x) = 
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— cosx — A cos 2x, hence L = 2n. With the choice A > 0, the direction of the particle current 
(X) turns out to be opposite to that of symmetry breaking in the SPDF P st (x): it is this 
effect which leads in turn to such oddities as negative zero-bias conductance and anomalous 
hysteresis Q]. The interaction force K(x — y) = ~K(y — x) between oscillators is a periodic 
function of x — y (also with period L) and in Refs. [f|,|5| is chosen as K(x — y) = K sm(x — y), 
with Kq > 0. We shall fix T = 2.0 and A = 0.15 as in Refs. 0,0], so the important parameters 
in the model are Kq and Q. The model just set up can be visualized (at least for A — > 0) 
as a set of overdamped pendula (only their phases matter, not their locations) interacting 
with one another through a force proportional to the sin of their phase difference (this force 
is always attractive in the reduced interval — 7r < x — y < tt). 

Mean-field approximation (MFA) With the above choice, the interparticle interac- 
tion term in Eq. (||) can be cast in the form 
1 N 

- ]T K(Xt - X 3 ) = K [d(t) sinX t - Si(t) cosX t ] . (3) 

7V 3=1 

For IV -> oo we may approximate Eq. (|3|) a la Curie- Weiss, replacing Ci(t) = 
iV -1 J2j cos Xj(t) by C m = (cosxj) and Si(t) = A^ _1 J2j sin Xj(t) by S m = (smXj), to be 
determined as usual by self-consistency. This decouples the system of stochastic differential 
equations (SDE) in Eq. ([I]), which reduces to essentially one SDE for the single stochastic 
process X(t). If X(t) happens to be Markovian, then it is a straightforward matter to write 
up an associated Fokker-Planck equation (FPE) whose stationary solution is the SPDF, 
from which all the transport properties can be readily obtained. As discussed in Ref . , the 
Gaussian character of the rji(t) allows in the r = case to consider them as being coupled 
through an effective function S(Xi) = V2[T + Q(W') 2 ] 1 / 2 . Hence in the white-noise case, 
the SDE obtained after performing the MFA is 

X = R(X) + S(X) V (t), (4) 

with R{x) = —V'(x) + F — A [C m sinx — 5 m cosx]. The associated FPE is 

d t P(x,t) = d x {-[R(x) + ±S(x) S'(x)]P(x,t)} + ^d xx [S 2 (x)P(x,t)] (5) 
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and its normalised stationary solution with periodic boundary conditions and current density 
J + is H 

p -<f>{x) lf( T ) 

where <j,(x) = -2 /* dy [ii(v)/S 2 W], H(x) = Jf" 1 * Sfe)" 1 exp [0(y)], and JV = 
J-i/2 dx P st (x). The positivity of S'(x) and the exponentials implies that of H(x) and hence 
that of P st (x) and AT, as it should be. 

The particle current: The appearance of a ratchet effect amounts to the existence of 
a nonvanishing drift term (X) in the stationary state, in the absence of any forcing (F — 0); 
in other words, the pendula become rotators in an average sense. As it was shown in Ref. 
||, the cause of this spontaneous particle current is the noise-induced asymmetry in P st (x). 
Being the current density 

J = [1 - e * L) \/2N, (7) 

the sign of J is that of 1 — e^ L \ The "holonomy" condition e^ L ' = 1 implies J = and 
H(x) = const. As shown in Ref. [HI it is 



(X) = [ L/2 dx \r(x) + ls(x)S'(x)] P st (x, C m , S m ), (8) 

J-L/2 L I 



with the result 



f i _ e <t>(L) ) 

{x)=JL -\-^\ L ' (9) 

hence (X) has the sign of J and can be also regarded as an order parameter. 

Equation (^) is a self-consistency relation since both M and <fi{L) carry information on the 
shape of P st (x) (in the latter case through C m and S m ). A nonzero J is always associated 
with a symmetry breakdown in P st (x) (namely, P st (—x) ^ P st (x)). This may be either 
spontaneous (our main concern here) or induced by a nonzero F. 

Effective Markovian approximation for coloured noise: The results in Eqs. (^) 
and (H) rely on the fact that we have been able to write the FPE Eq. (f|). When the 



i]i(t) in Eq. ([J) are coloured, the process X(t) in Eq. (f|) is in principle not Markovian 
and for non-Markovian processes, a FPE can at most result from some (non-systematic) 
approximation, like the truncation of some short correlation-time expansion. Fortunately, 
a consistent Markovian approximation (called "unified coloured-noise approximation" or 
UCNA) can be performed under certain conditions. By resorting to it one can obtain 
expressions for R(x) and S(x) in Eq. (gD which account for the effect of r. Their functional 
forms will be published elsewhere f|. 

The self-consistency equations: The stationary probability distribution P st (x) also 
depends on S m and C m , since R(x) contains these parameters. Their values arise from re- 
quiring self-consistency, which amounts to solving the following system of nonlinear integral 
equations: 

[L/2 

F cm = C m , with F cm = (cos x) = / dx cosx P st (x, C m , S m ), (10) 

J -L/2 
?L/2 

F sm = S m , with F sm = (sin a;) = / dx sinx P st (x,C m , S m ). 

J -L/2 

These equations give C m and S m for each set of the parameters (Q, K ) that define the state 
of the system, assumed T, A and F fixed. 

For F = 0, the choice S m = makes R(x) an odd function of x; this in turn makes 4>(x) 
even, and then the periodicity of P st (x) in Eq. (§) [in the formP s *(— x) = P st (—x — L)] 
implies that the stationary probability distribution is also an even function of x. So the 
problem of self-consistency reduces to the numerical search of solutions to Eq. (|T0"D, with 
S m = 0. Although plausibility arguments, detailed in Ref. 0, allow to have an intuition 
on the existence of some solutions to this integral equation (and their stability) in this 
symmetric case, the stability of the true solutions must be explicitly checked. Since cosx in 
Eq. ( |T0"D is an even function of x, it suffices to use the Curie- Weiss one-parameter criterion, 
namely to check whether the slope at S m of the integral in Eq. ( p"T|) is less or greater than 
one. As a complementary check, a small- a; expansion of <p(x) confirms that P st (x) is 
indeed Gaussian at x = 0. For small F^O, P st {x) gets multiplied (in this approximation) 
by exp[Fx/T] (= 1 + Fx/T) which leads to a nonzero value of S m = kF, with k > 0. By 
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(11) 



the mechanism discussed in Ref. ||, for large enough Q it is <f>(L) > and by Eq. ([?]) it is 
J < 0. This effect manifests itself in a negative zero-bias conductance since according to Eq. 
(|), (X) = LJ. 

As a consequence, for F = there are always one or more solutions to Eqs. ( |T0| ) and 
(|TT|) with S m = and one of these is the stable one in the "disordered" phase. As argued in 
Ref. |4j] , for N — > oo a noise- induced nonequilibrium phase transition takes place generically 
towards an "ordered" phase where P st (—x) ^ P st (x). In the present scheme this asymmetry 
should be evidenced by the fact that the solution with S m = becomes unstable in favor 
of two other solutions such that P 2 st ( x ) = A s *(~ x )-> characterised by nonzero values ilS^. 
This fact confers on S m the rank of an order parameter. 

The phase boundary: Since sinx is an antisymmetric function, Eq. (|TT| ) results imprac- 
tical for the task of finding the curve that separates the ordered phase from the disordered 
one, given that on that curve S m is still zero. For that goal (exclusively) we solve, instead 
of Eqs. fllUD and (|TT|), the following system: 

r£ / 2 , ~ „s ~ f L/2 . . dP st 



dx cosxP st (x, C m , 0) = C m , / dx sinx 

-L/2 J-L/2 Os 



= 1. (12) 

s m =o 



III. NUMERICAL RESULTS 



Figure |I] displays (in the same scale as in Ref. 0) the phase diagram obtained by solving 
Eqs. (|12D by the Newton- Raphson method. In the region enclosed by the thick lines ( "ordered 
region") the stable solution to Eqs. ( fLOf) and (|TTD has S m ^ 0. For r not too large this noise- 



induced phase transition is reentrant as a function of Q, for K = const (a fact already known 
for r = The novelty is that for any r^O, the phase transition is also reentrant as 

a function of Kq for Q = const (a feature found in Ref. || for a different system). 

The multiplicity of mean-field solutions in the ordered region, together with the fact that 
some of them may suddenly disappear as either K or Q are varied (a fact that is closely 
related to the occurrence of anomalous hysteresis) hinder picking out the right solution in 



7 



this region. A more systematic characterization of the aforementioned multiple solutions is 
achieved when the branch to which they belong is traced from its corresponding "homoge- 
neous" (S m = 0) solution. Accordingly, the thin lines in Fig. [I] separate two sectors within 
the ordered region with regard to the homogeneous solutions. Below them ("noise-driven 
regime" or NDR) there is a single solution with S m = and C m < (as already suggested, 
in this regime a solution with C m < can be stable since it corresponds to shaking violently 
the pendula). Above them ("interaction-driven regime" or IDR) there are three solutions: 
two of them have opposite signs and (for Kq/Q large enough) \C m \ ~ 0.9; the remaining one 
has C m ~ 0. Note that this line presents a cusp whose meaning was discussed in Ref. ||, 
in relation with the character of the hysteresis loop. We have studied the shape of P st (x) 
and the behaviour of (X) as a function of F for different locations in this (Q, K ) diagram. 
The square in Fig. [I] indicates a position inside the ordered zone for which the shape of the 
SPDF and the hysteresis cycle are followed as functions of r. This point lies in the NDR for 
t = and in the IDR for r / (marginally so for r = .1). 

Figure |2| shows (for the true solution, namely the stable S m ^ one) the evolution of 
function of r, for the state indicated by the square in Fig. [I] (Q — 10, K = 10.2). 
The SPDF is always an asymmetric function of x, indicating a spontaneous breakdown of 
parity (since V and W remain symmetric): the system has to choose between two possible 
asymmetric solutions, of which just one is shown. As discussed in Ref. ||, P st (x) is bimodal 
in the NDR. As r increases, it becomes unimodal (this is similar to what one achieves by 
decreasing Q at r = 0). 

Figures |3|(a) to ^|(c) present, for the state indicated with a square in Fig. [I], a sequence 
of (X) vs F plots obtained as r increases. All the solutions to Eqs. (|10"D and (|TT|), except 
the one belonging to the branch starting at C m ~ for S m = 0, have been included in these 
figures. The sequence is analogous to the one depicted in Fig. 9 of Ref. |J. In both cases 
we see a crossing from the NDR to the IDR (increasing r at fixed Kq has similar effects as 
increasing K at r = 0, for fixed Q). 
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IV. CONCLUSIONS 



From the analysis of the foregoing results, we draw the following conclusions: 

1. The increase of r tends to destroy order, and between r = 0.1 and r = 0.5 there exists 
a new reentrance with respect to K Q . In this range of r, the existence region of the 
ordered state is strongly shrunk. The qualitative similarity with the result arrived at 
in Ref. |J — in spite of the fact that both systems are different — shows the robustness 
of this result. 

2. Regarding the (X) vs F plots one sees that as r increases the hysteresis cycle becomes 
more complex, and the range of values of (X) corresponding to F ^ is severely 
limited. 

Although all of our results stem from a MFA, we see that this approximation is able to reveal 
the richness of the phase diagram of this model. Moreover, the mean-field results coincide 
with the numerical simulations we are undertaking and that will be published elsewhere [|J. 
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FIGURES 

FIG. 1. Phase diagram of the model for T = 2.0, A = 0.15, and F = 0.0. Full lines: r = 0.0; 
dashed lines: r = 0.1; dotted lines: r = 0.3. For each value of r, the ordered region lies above 
and to the right of the corresponding thick line. Above the thin lines there may exist several 
solutions when S m / 0, whereas below them there may exist at most one. The square corresponds 
to Q = 10.0, K = 10.2. 

FIG. 2. Shape of P st (x) at the point marked with a square in Fig. 1 (Q = 10.0, K = 10.2), 
for t =0.0 (full line), 0.1 (dashed line), and 0.3 (dotted line). The other parameters as in Fig. 1. 

FIG. 3. The order parameter V m = (X) (particle current) as a function of F for Q = 10.0 
and Kq = 10.2 (the square in Fig. 1). The sequence illustrates the change in the character of the 
hysteresis cycle as r varies: (a) r =0.0, (b) r =0.1, (c) r =0.3. Remaining parameters as in Fig. 1. 
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